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ABSTRACT 

' Context. The study of the snow line is an important topic in several domains of astrophysics, and 

particularly for the evolution of proto-stellar environments and the formation of planets. 
Aims. The formation of the first layer of ice on carbon grains requires low temperatures com- 
pared to the temperature of evaporation (T > 100 K). This asymmetry generates a zone in which 
bare and icy dust grains coexist. Our aim is to derive the proportion of bare grains around the 
theoretical snow line position for a typical low-mass protostellar disk and a massive protostar, 
and estimate the size of this mixing zone compared to the dust envelope size. 
' Methods. We use Monte-Carlo simulations to describe the formation time scales of ice mantles 

CO 



on bare grains in protostellar disks and massive protostars environments. Then we analytically 



describe these two systems in terms of grain populations subject to infall and turbulence, and 
assume steady-state. Numerical applications use standard numbers obtained by previous obser- 
vations or modelling of the astrophysical objects studied. 

Results. Our results show that there is an extended region beyond the snow line where icy and 
bare grains can coexist, in both proto-planetary disks and massive protostars. This zone is not 
negligible compared to the total size of the objects: on the order of 0.4 AU for proto-planetary 
disks and 5400 AU for high-mass protostars. Times to reach the steady-state are respectively es- 
timated from 10 2 to 10 5 yr. 

Conclusions. The presence of a zone, a so-called snow border, in which bare and icy grains co- 
exist can have a major impact on our knowledge of protostellar environments. From a theoretical 
point of view, the progression of icy grains to bare grains as the temperature increases, could be 
a realistic way to model hot cores and hot corinos. Also, in this zone, the formation of planetes- 
imals will require the coagulation of bare and icy grains. Observationally, this zone allows high 
abundances of gas phase species at large scales, for massive protostars particularly, even at low 
temperatures (down to 50 K). This could be a critical point for the analysis of upcoming water 
observations by the Herschel Space Observatory. 
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1. Introduction 

Inside the multiple tools that constrain the early phases of the formation of stars and plan- 



Universe, with «Hio/"Ht ~ 3 x 10 4 in the hot parts of molecular clouds ( 


Cernicharo et al. , 


1990: 


van Dishoeck & Helmich, 


1996 


Harwit et al. 


1998 


Boonman et al. 


2003 


) where it desorbs from 



the dust grain surface. Several observations of these hot parts have s hown that water plays a major 
role in the formation of complex molecules (e.g lBelloche et alj, l2009). This is due to the high dipole 



moment of this molecule that makes it active in both gas phase and surface chemistry. Water also 
acts as a host for complex molecules, keeping them in the icy mantle that covers d ust grains. In ad 
dition, water plays a role in the energetic balance of protostellar environments (e.g. 



Ceccarelli et al 



1996; 



Dotv&Neufeld . 



1997). The formation of planetesimals originates from the coagulation of 
icy dust grains: the formation of planets around stellar systems is triggered by the very first phases 
wher e dust grains stick to each other to form the bricks of rocky planets or the bulges of gaseous 



ones (iBrauer et al 



2008). Going a step further, water is the source of planetary oceans and is a 
necessary element for the appearance of life as known on the Earth. 

Previous results on the sublimation of dust ices from solid surfaces were made in laboratories, 



showing that it occurs on short time-scales from 100 K up to 160 K (iSack & Baragiola . 



Brown & Bolina, 



1993 



2Q07), and that the chemical products enclosed in the ice mantle are released si- 



multaneously. An astron omical example is given by complex molecules detected in h ot corinos and 



hot molecular cores (see 



Cazaux et al 



2003 



Ceccarelli et al 



2007 



Beuther 



2007 



for a review). 



Whereas the sublimation of ices is roughly known, its condensation or direct formation on grains 
is poorly understood. Ice mantles originate from accretion of gas phase water onto dust grains 



(iTielens & Hagen . 



1982: 



Cupnen & Herbs 



2007), or from direct formation of water via O and H 



chemistry on grain surfaces (ICazaux et al .1120 101) . Experiments and theory show that on graphite 
surfaces, water molecules a rrange themselves in clusters that have b inding energies that increases 



with the size of the cluster (ILin et al 



2005 



Brown & Bolina, 



2007). Once these clusters are cre- 



ated, water molecules are accreting onto themselves to form archipels from which the first layer 
of ice is made. This route is triggered by the condition that two molecules can meet each other on 
the surface to initiate the cluster. This p rocess is limited if the te mperature is too high, even for 
very high densities. In the case treated bv lCuppen & Herbstl (12007b . T ~ 10 K and ne, = 10 4 cm" 3 , 
appearance of the first monolayerof ice takes at least 10 3 years. This might be significant compared 
to a disk or a massive protostar evolution time-scale (from 10 5 to 10 7 yr), and where density and 
temperature conditions are different. 

At first, sublimation and condensation of water on dust particles has been considered to occur at 
the so-called snow line position, i.e. where the temperature in the protostella r environment reaches 
temperatures around 100 K (e.g. lSasselov & LecarluOOOtlLecar et ajj|2006j) . New results obtained 



on the formation of water ice show that the condensation of water is not optimal near this position, 
even impossible. Bare grains need to go back to a colder e nvironment to c reate the very first ice 
layer and which allows the ice mant le to accrete very fa st (Papoular, 



Padoan et al 



2002; 



MacLow 



2005). Protostellar environ- 



20041) . leading to the idea of an 



ments are known to be turbulent (e.g. 
extended zone where bare and icy dust grains can live together, creating the so-called snow border. 
Its extension depends on the material transport time-scales compared to the condensation times. In 
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this letter, we develop this idea for proto-stellar environments. We first investigate the time scales 
required for the formation of the first layer of ice. Then we derive standard equations that gives 
the mixing ratio ({r) between bare grains and the total amount of dust grains that can be expected 
around the classic snow line limit. We apply our model to a low-mass and a high-mass protostellar 
environment and conclude on the existence of a snow border. At last, we discuss how the snow 
border could affect present and future observations or modelling of proto-stellar environments, and 
the formation of planetesimals. 



2. Water accretion modeling 



We use a step by step Monte Carlo simulation to describe the formation of water cluster on surfaces. 
The dust grains are divided with square lattices of 100 x 100 adsorption sites. Each site of the 
grid corresponds to a physisorbed site where water molecules can be bound with a determined 
binding energy. We consider that water molecules accrete from gas phase, and that the rate of 
accretion depends of the density and the velocity of the water molecules in the gas phase. This rate 
of accretion can be written as: 



Ra 



nH,oVmooS 



(1) 



where n^o an d "HiO are the densities and thermal velocities of the water species, <x the cross 
section of the dust and S is the sticking coefficient of the species with the dust. We consider 5=1, 
meaning that when a water species arrive on a point of the grid, it becomes physisorbed. The water 
molecule that arrives on a location of the surface is adsorbed with an energy E. d that depends on the 
numbers of water molecules in neighbouring sites (see Table[TJ. Once the water is adsorbed on the 
surface, it can move randomly with a diffusion rate that can be written as 



fldiff = vexp 



[-¥) 



(2) 



where E&g = 2/3 E a is the energy of the barrier between two sites (from lCazaux & Tielens , 



2004) 



T the temperature of the surface and v the oscillation factor of the water molecules (typically 
10 12 s _1 for physisorbed species). Adsorbed molecules can also evaporates with a rate /? e vap iden- 
tical to Raff except that E&g is replaced E- d . The molecules on the surface can execute a random 
walk on the surface that generates the formation of water clusters. As two water molecules are in 
neighboring sites, a oxygen bound interaction link the two species, making the binding energy of 
the dimer higher than the binding energies of individual molecules. Then, other water molecules 
can join and increase the size of the cluster, and, as a result, in crease the total binding energy 
( Gonzalez et alj|2007 ). We used the binding energies calculated in Lin et al. ( 2005 ). as reported in 



TableQ] A water cluster grows along the surface until it reac hes six members. Whe n the size of the 
cluster is larger than six molecules, the cluster grows in 3D ([Gonzalez et al.l 120071) . with a binding 
energy dependin g only on the number of water molecules which are close to the graphite surface 



dLin et al 



20051) . Therefore, in our model, we limit the number of water that compose a cluster to 



6 water molecules. The formation of cluster is therefore a competition between the evaporation of 
individual water molecules and the encounter of two molecul es to initiate the cluster . In this work, 
we only consider carbon surface. However, recent sutudies by|j 



Goumans et al 



(2009) show that the 



adsorption of water on silicate grains has very high binding energy (10000 K). Therefore, for these 
type of grains, we don't expect an assymetry between ice formation and evaporation. 
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Table 1. Adsorption energ y (in kelvins) of water clusters as a function of the number of molecules 



present in the cluster, from 

n=l n=2 n=3 



Lin et al.B2005 



n=4 



n=5 



n=6 



E„ [K] 1440 2400 3400 4600 5400 6400 
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n b (r) 



K" 

n b (r+dr/2) 
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r-dr/2 r+dr/2 



Fig. 1. System description used for the analytic description of the dust grain movements. We con- 
sider very small elements of width Ar in which density of bare grains are constant with the value 
rib(r). Between a small time range Af, due to turbulences, diffusion of bare grains («i(r)) is possible, 
triggered by K^r) fluxes (see Sect. 3). 



3. The mixing ratio 

In order to analytically describe the system, we consider a small zone around a radial position r (see 
Fig [TJ with an extension Ar. Inside this zone, all physical parameters are assumed to be constant. 
Going to infinitesimal variations in time and space then allow us to derive the equations describing 
the system evolution. We follow the bare grain population, which evolves through the differences 
between the fluxes of grains moving away from the central star K + and the fluxes of grains going in 
the direction of the central star Kr. These fluxes are obtained from the description of the velocity 
fields inside the considered protostellar environments. To study low-mass proto-planetary disks and 
massive protostar environments, we consider cylindrical and spherical descriptions of the system. 
For both of them, we derive the mixing ratio £(r) between bare and total density of grains, as 
the system reaches steady state. Thus this mixing ratio describes the snow border itself. The first 
subsection gives the analytic description of the systems and the velocity fields acting on them, then 
the next subsection discuss the time scales and feasibility of reaching steady state. 



3.1. Velocity fields 

We consider three main components acting on the velocity field: the thermal velocity Vg, the turbu- 
lent velocity u t and the accretion velocity u acc . The two first components are assumed to be isotropic 
and can be averaged geometrically to form a composite velocity u T : 



1>T 



2k B T 2 



1/2 



(3) 




Fig. 2. Overview of the snow border in a proto-stellar environment. At a certain distance from 
the central star, there is a zone where bare and icy grains coexist. This zone is located between 
two regions: maximum temperature region at which ices evaporate (here 160K) and a minimum 
temperature region at which the first layer of ice can form (here 60K). Turbulences allow grain and 
gas to mix. 

where T is the dust temperature, k% is the Boltzmann constant, m is the average mass of a dust 
grain and v t can be determined with molecular line observations. The corresponding velocity field 
can thus be assumed to follow a Maxwell-Boltzmann distribution of the form: 

f{v) = %v 2 e- aV (4) 

71 1 

where a = 1 /uj. In our case, we are interested in the radial distribution of the velocities which 
allow the transport of grains to higher or lower temperatures. In proto-stellar environments, where 
the accretion of materials onto the central protostar is done radially through disks or shells, the 
radial field f r (u) follows the form: 



2,2 



Mv r ) = -^e-™' (5) 

This velocity distribution is shifted by the velocity field induced by the accretion. Knowing the 
accretion rate M of the protostellar object permits to deduce v^c- 
M 

UaccOO = —. . r , cylindrical, (6) 

47Trfim p n(r)H(r) 

M 

= —. — 5 tt> spherical, (7) 

4nr~jim p n{r) 

where is the average mass of a gas particle, m p the proton mass, n{r) the IT density and H{r) 
the height scale of the disk at the distance r. We note that vj also depends on r since it contains a 
thermal component varying with T(r). The new radial velocity distribution f*(v r ) is then obtained 
by shifting the initial distribution: 

f*(Vr) = fr(Vr ~ "ace) (8) 

This result allows us to define an average velocity for dust grains travelling to lower radius, denoted 
(v~), and another one for the opposite trend, called (v + ): 

<if>(r) = f vf?(v r )du (9) 

%J —CO 

r*+oo 

(u + )(r) = vf;(v r )dv (10) 

Jo 
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giving 

<OW = ^[1 -Erf(au acc )] - (11) 
<» + > W = [1 + Erf ( fl! ; acc )] + ^— ^ (12) 

where Erf(x) is the standard error function. Now that the mean velocities of inward or outward mo- 
tions are known, we derive the fractions of dust grains A^ ± (r) having positive or negative velocities. 
They are given by: 

N-(r) = f f;(v r )dv (13) 

tj CO 

f;( Vr )dv (14) 



hence 

AT(r) = i[l-Erf(a(rK cc (r))] (15) 

A» = ±[l+Erf(a(r) Uacc (r))]. (16) 

These numbers allow to derive flux densities of dust grains at a certain position r in the model 
that we define. These fux densities depend on the geometry considered: 

^(r) = +27rrtf(r)<u ± )jV ± (r), cylindrical, (17) 
= +7rr 2 <i; ± )jV ± (r), spherical, (18) 

where the sign + in front of the expressions are introduced to make the value of flux densities 
always positive. Indeed, from its definition, (tr) (r) has negative values. 

3.2. Radial dependence: £ (r) 

We consider the variation in bare grain density «/,(r) inside a shell of width Ar around the position 
r (see Fig. ??). Bare grains move from r with flux densities ^T ± (r) or disappear into icy form with 
a speed k c (r) (in grains. s _1 ). Determining this last parameter in then possible through the modeling 
we described in previous section, and empirical formulae will be given for each case studied in the 
next section. Indeed, k c depends on r as it is linked to the local conditions, i.e. n(r) and T(r), and 
can present huge variations depending on the proto-stellar environment studied. 

Variations in the total number of bare grains Ni,{r) between r - Ar/2 and r + Ar/2 are then 
described by: 

N b (r, t + At) = N b (r, t) - N bomn + N in - N out (19) 

where Nf Iozen is the amount of bare dust grains becoming frozen, i.e. having formed 10 % of the 
first ice layer, per time unit. This number is related to k c by: 

AW„ = k c (r)N b (r,t)At (20) 

but can also be related to the local density introducing the volume element depending on the model 
geometry: 

AW en = K c (r)n b (r,t)ArAt (21) 
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hence 

K c (r) = AnrH(r)k c (r), cylindrical, (22) 
= 2jrr 2 k c (r), spherical. (23) 

Other terms in Eq. [19] N- m and N out , refer to the incoming and outcoming numbers of grains due to 
velocity fields described previously. These fields are inducing fluxes that depend on the geometry 
considered: 

K^r) = ±2nrH(r) <i; ± > N*(r), cylindrical, (24) 
K^r) = +7rr 2 <y ± )/V ± (r), spherical, (25) 

giving the following expressions: 

+ K-(r+^)n b lr+^-)At (26) 



2 ) \ 2 

N out = K + (r)n b (r)At + K-(r)n b (r)At. (27) 

Going to the limits when Af — > and Ar — » yields the differential equation 

— (Jn b ) - n b K c = (28) 
ar 

with 

/(r) = * + (r) - AT(r) (29) 

which can be solved numerically with the boundary condition that all the ice mantles are totally 
sublimated beyond 160 K (n b (r s ) = n(r s )). This temperature corresponds to the de sorption of the 



most strongly bound ice st ructures that can be cre ated on dust (i.e. thick water ice. ICollings et al 



2004, on graphite surface, 



Brown & Bolina 



2007). The snow border extend can then be derived 



through the ratio ( (r) given by: 

^) - (30) 

We then assume in our study that the snow border extension AR corresponds to the distance between 
the ice mantles sublimation radius and the radius where f(r) = 0.5. 



3.3. Time estimate to reach steady-state : t 

The snow border extend refers to the steady-state of the system. In this sense, one must consider 
the time to reach this equilibrium. In our study, an estimate of this time is critical, as proto-stellar 
environments have a limited lifetime linked to the accretion rate and the final star mass which is 
formed. As the system must be stable on a given period to fit our assumptions on the density and 
temperature distributions, the snow border idea makes sense only if its growth is fast enough. 

To roughly estimate the time t needed to reach the steady-state, we focus on the mean velocity 
of grains going to the external part of the proto-stellar environments, i.e. (v + ). An estimate of t is 
then given by: 
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As this value depends critically on the diffusion length scale in real systems, this estimate assumes 
that turbulences are intermediate scale motions compared to the size of the proto-stellar environ- 
ment. This makes this first estimate of r a lower limit to reach the steady-state, as the diffusion 
speed can critically increase if the turbulence typical length scale is very small. In our case, disk- 
like environments are the most difficult to treat, and t estimates we obtain can be realistic restricting 
the model to the mid-plane, where turbulences are near to the thermal width, and gas is coupled 
with grains. For high-mass protostars, both study of molecular clouds dynamics and observation 
gives a correlation leng t h between 0.05 pc and 0.1 pc, i.e. the typical s ize of a massive dense core 



(Miesch&Ballv 



1994 



Goodman et al 



1998 



Ossenkopf et all 120011) . We can then be relatively 



confident with this first estimation, which is also coherent with the first-guess analytic approach we 
have in this study. 



4. Astrophysical applications 

We first apply the model described previously to two particular proto-stellar environment where 
the snow border may have an important extend, hence impact. We choose a typical proto-planetary 
disk and a massive dense core, in which low-mass and high-mass star formation occur. This goal 
of this first step is to disantangle the fundamental parameters that trigger the snow border extend 
and growth, before going to a general test on the model covering a grid of these parameters. 



4.1. Proto-planetary Disks 

The snow line position is very important to constrain the formation of planetesimals in low-mass 
star forming regions, as current models treating this topic take into account bare-bare or icy-icy 
grain collisions only. During the formation of low-mass stars, accretion of materials on the cen- 
tral protostar is done through a rotating gaseous and dusty disk, which is the precursor of the 
planetary system. As a first test case for our modeling, we u se for our application a model of 



Woitke et al. 



(2009). Turbulence, observed with 



a proto-planetary disk around a T-Tauri star by 
molecular line emi ssions towards a d isk surrounding the Herbig Ae star HD 34282, is of the order 
of 0.1 km s _1 (e.g. 



Pietu et al 



20031) . This value, has been estimated for a Herbig AE disk that is 



more massive and somewhat hotter than the TTauri disk considered in this study. Since such turbu- 
lent bandwidths have not been detected towards TTtauri disk, we consider 0.1 kms -1 as an upper 
limit. We also consider a mas s accretion rate in this object of 10~ 6 M Q yr~', which is also an upper 
limit ( Hueso & Guillol2005 ). According to the turbulence limitation we discussed in the previous 
section, we restrict the snow border extend and growth in the mid-plane of the accretion disc. This 
also implies a very low UV field for the modeling of the water accretion on carboneous dust grains. 
In the mid-plane, H2 density is known to follow a power-law of the form: 

\p 



n{r) 



'0 



(32) 



where «p is the density at the reference radius ro and p a power-law index, equal to -2.45 in the 
model by 



Woitke et al 



(120091) . The disk height scale H{r) is proportionnal to r 2 and the temperature 
T(r) is proportionnal to r~ 025 , both of them following again a power-law similar to the one indicated 
for the density. Values of reference at ro = 1 AU are «o = 10 14 cirT 3 , Tq — 130 K and Ho = 0.1 AU. 
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Table 2. Results of the water accretion modeling on carbenous grains in the mid-plane of a proto- 
planetary disk (left) and in high-mass proto-stellar object (right). 



Conversion to dust grain density assumes a standard MRN distribution (IMathis et al.lll977l) with 
gas-to-dust ratio of 100. 

Five density-temperature couples corresponding to the proto-planetary-disk are used to derive 
an empirical formula giving the condensation speed k c , from the results obtained by the grain sur- 
face modeling. The water abundance assumed, relative to H2, is 10 -4 according to atomic universal 
abundances. Ice mantle growth rates obtained are given in Fig.[3j and the corresponding condensa- 
tion speeds are given in details in Table [2] Resulting formula for the proto-stellar disk case in the 
mid-plane is: 

k c (r) a 6.65 x 10"V - 12 * r(r) [s -1 ]. (33) 

In this case, the limit where T > 160 K is r eva p = 0.43 AU. The snow border extension is shown 
in Fig. [4] The result shows an extended region where the mixing ratio is significant. The domain 
where f > 0.5 is spread over Ar ~ 0.2 AU, that corresponds to T ~ 135 K. A first conclusion is that 
the snow border does not reach its maximally possible extent, because diffusion of bare grains is 
opposed by the accretion velocity field and re-formation of ice mantles. Evaluation of the diffusion 
time t ~ Ar/ (v + ) (r evap ) gives the steady-state value r ~ 10 2 yr, i.e. shorter than the disk lifetime. 

4.2. High-Mass Protostars 

The snow border is particularly interesting to follow in the case of massive protostars, where strong 
turbulence is observed (v t > 1 kms~') in addition to high accretion rates (M > 10~ 5 M Q yr -1 ). 
Furthermore, due to their high luminosities (L > 10 4 L ), the initial snow line position in such 
objects is shifted to higher distances. At this location, gas density drops by orders of magnitude 
compared to low-mass protostars (10 5 to 10 6 times lower) which should have a strong effect on the 
formation of the first layer of ice. 

In this study we consider a turbulent velocity field value of ~ 1 km s , as shown by molecular 
line width measurements and modeling. We also include a mass accretion rate of lxlO -3 M yr -1 
as required for the formation of high-mass stars, and consider the spherically symmetric expres- 
sions f or and k c . The earl y stages of high-mass protostars have been modeled considering a ID 



shape ( Marseille et al. 2008), which allow us to derive the gas density and temperature as function 
of the distance from the central star. The gas density distribution used follows a power-law similar 
to Eq. [32] with ro = 21 AU, «o = 2.8xl0 9 cirT 3 and p - —1.2, giving a temperature distribution 





10 -9 I , 1 , , , 1 , , , 1 , , , 1 — 

60 80 100 120 

T [K] 



Fig. 3. Plot of the bare-to-icy grain transformation rate k c as function of the temperature in the 
proto-planetary disk considered. 



following the same behavior, with To = 1500 K and p - -0.59, accord ing to the modeling of the 



IRAS 18151-1208(MM1) massive dense core in 



Marseille et al 



(2008). 



As for the case of proto-planetary disk, we define several density-temperature couples that 
corresponds to high mass protostellar objects. For each of these couples (see Table|2), we calculate 
the time to form the first layer of ice using our grain surface modeling. With these calculations, 
we can derive an empirical formula giving the condensation speed k c . In a similar way than for 
proto-planetary disks, we consider the water abundance of 10~ 4 relative to H2. Ice mantle growth 
rates obtained are given in Fig. [5] and the corresponding condensation speeds are given in details 
in Table[2] Resulting formula for high mass protostellar objects is: 



k c (r) - 2.97xl0- 8 e - a23xrw [s- 1 ]. 



(34) 
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Fig. 4. Plot of the mixing ratio £(r) as a function of the radius for: T-Tauri proto-planetary disk 
models (solid line, bottom labels), massive protostars (dashed line, upper label). 

The derived mixing ratio ({r) shows that the domain where ( > 0.5 is spread over ~ 5000 AU 
(see Fig. |4), reaching T ~ 55 K. However, the diffusion time scale is much longer in this case 
(t ~ 10 5 yr), indicating that the system may not have time to reach this steady state completely. In 
the case of high-mass protostars is then very extended because of the high turbulences and the low 
densities in the region where water freezes-out onto dust. 

5. Impact of the main parameters 

The results of the previous section have been obtained for what we assume to be "typical param- 
eters" in the considered proto-stellar environments. Although a clear difference between low-mass 
and massive protostars is present, the dependence of the snowborder caracteristics (extension and 
growth speed) on the main parameters should be explored. These main parameters, linked to the 
mixing effect which is mandatory for the snow border formation, are the stellar accretion rates and 
the turbulent velocity fields. 

Keeping the proto-planetary-disk and the massive protostar as basis for the following study, 
we now explore a larger range of value for v t and M. For proto-planetary-disks, the turbulent ve- 
locity field ranges between 0.01 and 0.1 kms~' and the stellar accretion rate between 10~ 8 and 
10~ 5 M Q yr -1 . Concerning high-mass protostellar objects, these values range between 1.0 and 
10 kms 1 for the turbulent velocity field and between 10~ 5 and 10~ 2 M yr 1 for the accretion 
rate. 

5.1. Impact on snow border growth 

The Fig. [6] (color scale) describes the influence of the stellar accretion rate and the turbulence 
velocity on the growth of the snow border in AU yr -1 . The snow border growth ranges between 
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Fig. 5. Plot of the bare-to-icy grain transformation rate k c as function of the temperature in the 
high-mass protostellar objects considered. 



and 0.005 AU yr~' for the proto-planetary model and between and 0.5 AU yr~' for the high-mass 
protostellar object model. In this sense, a null value means that the snow border formation will not 
occur. The snow border growth increases as stellar accretion rate decreases and turbulent velocity 
increases. On one hand, the stellar accretion prevents the diffusion of water molecules from the 
inner part of the objects to higherradii. On the other hand, turbulences support the diffusion towards 
the external parts. As a consequence, accretion rates lower than 10~ 5 M yr~' (resp. 10~ 2 M yr ) 
allow the appearance of the snow border in protoplanetary disks (resp. HMPO). At the same time, 
turbulences higher than 0.2 km s _1 (resp. 2 km s _1 ) are required for these environments. 
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5.2. Impact on steady-state time 

Fig. [6] (labelled white contour map) shows the time to reach the steady state equilibrium. This 
time increases strongly with accretion, and slightly with turbulent velocities: accretion slows down 
diffusion while turbulence improves it. The appearance of an extended snow border can only occur 
if its associated growth is important enough and the time to reach the steady-state is sufficient. 

In the case of proto-planetary disks, the maximal extension for the snow border is obtained 
after 100 years for turbulence of 0.1 kms 1 and an stellar accretion rate of 5xl0~ 6 M yr _1 , and 
reaches 0.2 AU (0.2 % of the total size). This time is reasonnable compared to standard evolution 
time scale of a protoplanetary disk, meaning that the snow border exists. 

In the case of high-mass protostellar objects, the maximal extension can be obtained after 
10 5 years, for a turbulent velocity field of 10 kms _I and an accretion rate of 5xl0~ 3 M yr _1 , 
and reaches 10000 AU (40 % of the total size). Above this time, HMPOs are known to enter a 
phase of medium ionization due to the central star starting its life on the main sequence. 

Nevertheless, the cases mentionned above are extreme compared to standard objects observed. 
More reasonnable value for the total size of the snow border should be around 0.04 AU for proto- 
planetary disks and 5000 AU for HMPOs. This values correspond to the lower left part of the graph 
in Fig. 

6. Discussion and conclusions 

The condensation of water ice under 100 K requires long time scales. This allows the coexistence 
of both icy and bare grains in a zone beyond the snow line, that can be called the snow border. This 
region can be significant compared to the size of the astrophysical objects considered in this work, 
particularly for massive dense cores where its extent can reach a fourth of the total system size. 
More precisely, the existence of the snow border is triggered by the stellar accretion rate and the 
turbulent velocity field in the protostellar object. Our study shows that a minimal turbulent level of 
~ 0.1 km s _1 in protoplanetary disks (resp. 1 km s in HMPOs) is required to allow fast diffusion 
of water molecules. It also indicates that stellar accretion should be lower than ~ 10~ 5 M yr~' 
in protoplanetary disks (~ 10~ 2 M yr 1 in HMPOs) to allow an efficient appearance of the snow 
border region, otherwise the water diffusion towards the object is prevented. 

The snow border seems critical for observations of massive dense cores, where its size (few 
thousand AU at few kpc, hence ~ 1") is not negligible. Indeed, to first order, the existence of such 
a region extends the zone in which molecules are present in the gas phase. Interferometers, and 
obviously the ALMA project in the coming years, should be able to reach the spatial resolution 
corresponding to the snow border extent. More precisely, the snow border allows a non-negligible 
abundance of water in the gas phase, even at cold temperatures (down to ~ 50 K). This could be 
a critical point for the analysis of the upcoming water observations in massive dense cores by the 
Herschel Spatial Observatory. However, we must stress that our study is based on carbon surface 
behaviours, whereas astronomical silicates react differently. Particularly, a fast condensation of 
water on silicate remains below 160 K. As a consequence, a sharp release of complex molecules 
over 160 K is still valid, whereas the snow border region will be characterised by a transfer of 
ice materials from carbon to silicate grains. In the future, this issue should be investigated, this 
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Fig. 6. Overlay of the snow border growth (color scale) and the time to reach the steady-state (white 
contours) as functions of stellar accretion rates and turbulent velocity field in protoplanetary discs 
(left) and high-mass protostellar objects (right). 



discrepancy between the two type of grains being possibly important in the physical and chemical 
evolution of protostellar environments. 

From a modelling point of view, this result opens new doors on several topics, and may have 
a non-negligible impact on the evolution of star forming regions. First, the snow border existence 
allows a smooth transition of molecular abundances in the gas phase and on the grain surfaces. This 
is new compared to the abundance jump assumed commonly in chemical modellings or observation 
analyses. A natural next step will be to couple gas and grain chemistry to better constrain these 
modifications. In addition, as the formation of the first layer of ice depends on the grain size, it 
is of prime importance to study this aspect in further detail. Concerning the modelling of strongly 
heated objects like of hot cores or hot corinos, the snow border may modify the grain temperatures, 
hence the results obtained. The coagulation of bare and large icy dust grains could also lead to 
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the creation of new type of proto-planetary objects, through collisions between these two types of 
grain s. Indeed, current modelling consider bare-b are or icy-icy grain in ter actions only (Ormel et al 



20071) . In addition, previous results obtained by 



Kretke & Lin (2007) or 



Brauer et al 



(2008) tend 



to show that formation of planetesimals is particularly enhanced in the neighbourhood of the snow 
line. Concerning the formation of planets, the fact that bare grains can exist at large distance from 
the central star could also play a role in the rocky vs gaseous planet formation issue, pushing its 
limit further away. 
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